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We study the superfluid state of atomic fermions in a tunable optical kagome lattice motivated 
f^*) ' by recent experiments. We show that imposed superflow induces spatial modulations in density 

£SJ ' and superfluid order parameter and leads to a charge modulated superfluid state analogous to a 

supersolid state. The spatial modulations in the superfluid emerge due to the geometric effect of the 
^ 1 kagome lattice that introduces anisotropy in hopping amplitudes of fermion pairs in the presence 

Qof superflow. We also study superflow instabilities and find that the critical current limited by the 
dynamical instability is quite enhanced due to the large density of states associated with the flat 
band. The charge modulated superfluid state can sustain high temperatures close to the transition 
temperature that is also enhanced due to the flat band, and is therefore realizable in experiments. 
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Geometric frustration is a central subject in modern condensed matter physics. Various ordered and liquid phases as 
well as multiferroic behavior can arise from geometric frustration. The kagome net is a well known example of lattice 
geometry that exhibits high degree of frustration. This lattice geometry is proposed to host various exotic phases such 
as quantum spin liquid, valence bond solid [H-[H , and fractional quantum Hall state [j| . The intriguing feature of the 
kagome lattice is the nondispersing flat band arising from geometric frustration. It enhances interaction effect and 
leads to ferromagnetic order [§[ as well as the destruction of Bose- Einstein condensation and the resulting supersolid 
state 0- Superconductivity (supefluidity) on the kagome lattice is a topic of great recent interest from theoretical 
side despite that few corresponding systems are known Q . The infinitely large density of states associated with 
' the flat band can strongly enhance superconductivity in the kagome lattice Geometric frustration also provides 
£h ■ highly nontrivial effect on electron correlations and thus can lead to the emergence of novel superconducting states 

m 

Ultracold atoms trapped within optical lattices offer an ideal system for explo ring various exotic phases and studying 



quantum phase transitions due to its remarkable controllability and cleanness [llj . Stability of superflow and critical 



velocity are of particular interest for Fermi and Bose superfluids in optical lattices since the pioneering experiments 



by Ketterle's group [12| . Various dynamical instabilities arising from competing orders in optical lattices and the 
q>^ ■ possibility of "flowing supersolid" state have been investigated 13l4l7j|. The experiments by Jo et al. [18( realized a 
£SJ ' tunable optical lattice in the kagome geometry by overlaying two triangular lattices with commensurate wavelengths. 

■ Introducing fcrmionic isotopes 6 Li and 40 K into optical kagome lattices, this system would provide a very important 
^vq | platform for studying superconductiviy and Cooper pairing in the kagome lattice. 

7-H . In this Letter, motivated by the recent experimental developments reported by Jo et al. [Isj] . we study the superflu- 
idity of atomic fermions in an optical kagome lattice within the attractive fermion Hubbard model. Imposing a nonzero 
superflow, we calculate superfluid order parameter as well as fermion density within the mean-field approximation. 
^ ' We further study superflow instabilities which limit the critical current by evaluating the free energy as a function 

J^j . of the flow momentum. Figure [1] shows the stability phase diagram of a s-wave superfluid state that summarizes the 
main results of the paper. We find a novel geometric effect of the kagome lattice that leads to a stable flow-induced 

■ charge modulated superfluid state where the superfluid order and charge density wave (CDW) order coexist (see the 



inset of Fig. [TJ) analogous to "flowing supersolid" proposed in Ref. |l3|. Furthermore, the critical current for this state 
is found to be quite enhanced at high fermion density due to the diverging density of states associated with the flat 
band. We also find that the instabilities at the critical current are dominated by the different mechanisms at the 
low and high fermion densities. The stability phase diagram thus shows a remarkable particle- hole asymmetry of the 
critical current. The unexpected charge modulations we uncover can be observed by em ploy ing the setup of the recent 



cold-atom experiments [l8| combined with the moving optical lattice technique 12 , hJ 2(j or dipole oscillations 21 
The charge modulations induced by superflow may be also observed around a quantized vortex. 

We consider two-component atomic fermions with attractive interactions in a deep optical kagome lattice. The 
system is well described by the attractive Hubbard model: H = — £$^(i ^ a {$i a Cjo + h.c.^ — U n^hii ~ M a n ia-> 
where c\ a is the creation operator of a fermion with spin cr(=y,4,) at site i, hi a = cj^c^, t(> 0) is a hopping amplitude 
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between nearest-neighbor sites, U(> 0) is the local Hubbard attraction, and the chemical potential /i tunes the 
fcrmion density. For U = 0, we can diagonalize the Hamiltonian in momentum space and obtain the three energy 



bands: £o = 2t and e±(k) = t y— 1 ± J 3 + 2 J2i=i 2 3 cos (k ' a i)J with ai = (1,0), a2 = (1/2, \/3/2), and a3 = a2 ai. 

The highest band £0 has a flat dispersion and the other two bands form Dirac cones at the K and K 1 points in the 
Brillouin zone. The density of states (DOS) D(e) of free fermions diverges at e = 2 due to the flat band and vanishes 
at e = —1 due to the Dirac points. The flat band is occupied if the fermion filling n is in the range of 4/3 < n < 2, 
while the Fermi energy is at the Dirac points when n = 2/3. 

We treat the on-site attractive interaction within the Hartree-Fock-Gor'kov (HFG) approximation. We in- 
troduce the superfluid order parameter in the presence of superfluid flow with flow momentum Q, = 
(U/M) Z)k( g -k+Q/2,i/,4.Ck+Q/2,i/,t)) and tlie average number of atoms per site, n v = (1/M) ^^(c^o-Ck^o-), where 
Ck.i/,cr is the Fourier component of c !( j, M is the number of unit cell, and v(=A, B, C) labels the sublattice as shown 
in the inset of Fig. [T] The fermion filling is given by n = (jia +ijb + nc)/3. Within the mean-field approximation, 
the Hamiltonian takes the form 

^HFG = E VI 'Q( k ) /l Q( k ) VI 'Q( k ) + £; Q' W 

k 

with 

*q( k ) = ( 6 k + ,A,t £ k_A,4 4+,B,t 6k -3,l 4+,c,t ^-,0,4.) 

and = -3M(/j, + Un/2) + MJ2A\^u\ 2 /U + Un 2 jA). Here, k± denotes ±k + Q/2. The explicit form of h Q (k) 
is presented in the Supplemental Material (22[. We obtain the excitation spectrum of Bogoliubov quasiparticles 
i?q jT (k) for the energy band r(= ±,0) as the eigenvalues of the matrix /iQ(k). The three gap equations for A v 
and the three number equations for n v are obtained from the unitary matrix for the Bogoliubov transformation that 



diagonalizes ^iQ(k) in the standard manner 17] . To evaluate A„ and /i, we solve these gap and number equations 




FIG. 1: (color online) Stability phase diagram of a flowing superfluid state in the kagome lattice. We set the flow momentum 
Q = (Q,0) and T = 0. Qd and Q C 2 are the critical flow momenta for dynamical instabilities with distinct mechanisms (see 
the text). The inset shows the stripe charge order induced by superflow with n\ = »b / n c (see the text). 
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FIG. 2: (color online) Density modulation h, order parameter modulation A, and averaged order parameter A as functions of 
flow momentum Q for (a) n = 0.4 and (b) n = 1.6 when the flow is in the V — K direction (Q = (Q, 0)). We set U/t = 3 and 
T — 0. The vertical dashed line in each panel marks the location of the critical flow momentum Q c ± or Q C 2 above which the 
system becomes dynamically unstable. 



self-consistently. This scheme is known to interpolate the weak-coupling BCS regime and strong coupling BEC regime 
at low temperatures [23j]. 

In the stationary ground state (Q = 0), the order parameters as well as the filling factors take the same value for 
all the sub-lattices by the symmetry, i.e., Aa = Ab = Ac and ua — ub = uq. If we impose nonzero superflow, we find 
that infinitcsimally small amount of superflow breaks this symmetry and leads to spatial modulations in density and 
order parameter. We plot n = nc — ua, A = Ac — Aa, and the averaged order parameter A = (Aa + Ab + Ac)/3 
as functions of Q for different fillings in Figs. [2] (a) and (b) setting the flow in the T — K direction (Q = (Q, 0)). The 
system involves a stripe order in the direction perpendicular to the flow with Aa = Ab ^ Ac and ua = =/= nc 
(see inset of Fig. [I} . This superfluid state with stripe charge order is analogous to a supersolid state in the sense 
that superfluid and CDW orders coexist. Such a flow-induced charge modulated state arises due to the geometric 
effect unique to the kagome lattice, while the supersolid state arises due to the spontaneous breaking of translational 
symmetry. 

To check the stability of the superfluid state with stripe charge order, we evaluate the quantity d Jq2^ , where 
N = 3M is the number of total lattice sites and F{Q) is the free energy of the system in the presence of superflow. Note 
that the superfluid density is proportional to d 2 F(Q)/dQ 2 \Q = Q. If this value is negative, the system is dynamically 

unstable against phase and density fluctuations [l~3|, EH- As shown in Figs. [2ja) and[2jb), d Jq£^ is positive for 
both the low and high filling until the critical flow momentum Q c \ or Q C 2- Consequently, the superfluid state with 
stripe charge order is dynamically stable until Q c \ or Q C 2- This is in sharp contrast with the "flowing supersolid" 
state in a square lattice discussed in Ref. [lj. This state with the superfluid and checkerboard CDW orders is found 
to be dynamically unstable with negative d Q2 for any nonzero Q. 

We find that the dynamical instabilities are caused by different mechanisms for low [Fig. [2] (a)] and high fillings 
[Fig. |5](b)]. The value of d in Fig.[2ja) changes discontinuously at the critical momentum Q c i, while the curve 
in Figj2jb) smoothly changes from positive to negative at Q C 2- We find that the sudden change of the curve at Q c i 
is associated with gapless quasiparticle excitations. With increasing flow, the energy gap to creating quasiparticles 
decreases due to a Dopplcr shift in the direction opposite to Q as shown in Fig. [3] (a). For large Q, the gap closes 
and the lowest quasiparticle energy band i?Q, T (k) becomes gapless. This precisely coincides with Q c \. Since the free 
energy at T — involves the contribution from spontaneously excited quasiparticles with negative energies above the 
critical flow Q > Q c i, the second derivative of the free energy ~ fg? changes discontinuously at Q = Q c \. On the 
other hand, in Fig.[2jb), the quasiparticle dispersion is still gapped at Q C 2- Thus, the instability at Q C 2 sets in before 
the closing of the single-particle excitation gap. The negative value of d therefore indicates the dynamical 



instability associated with collective phonon excitations rather than single-particle excitations [l3j, l2J-|26[ . At the 
onset of the instabilities at Q cl and Q C 2, the frequency of long- wavelength phonons becomes complex. As a result, 
the amplitude of collective phonon excitations grows exponentially and the superfluid state collapses. 

Figure [T] shows the critical momentum Q c as a function of filling factor n when the superflow is in the r — > K 
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FIG. 3: (color online) Quasiparticle energy band _Eq >T (k) (r = 0, ±) for different values of imposed flow Q at (a) n — 0.4 and 
(b) n = f.6 when the flow is in the T — K direction (Q = (Q,0)). We set U/t = 3 and T = 0. The critical momenta are 
Qd « 0.31tt for (a) and Q c2 ~ 1.04tt for (b). 



direction. It exhibits a remarkable particle- hole asymmetry reflecting the different features of at Qci an d 

Q C 2 discussed above. The critical momentum at high filling (n > 4/3) is significantly enhanced being limited by the 
onset of the dynamical instability at Q C 2- The robust superfluidity against imposed superflow is due to the flat band 
in the non-interacting band structure. The diverging DOS at the flat band enhances the superfluid order parameter as 
well as the single-particle energy gap [compare Figs.[3][a) andEJb)]. The large energy gap requires large Q for closing 
it by a Doppler shift and thus suppresses depairing instability. On the other hand, the small order parameter for low 
filling yields the small energy gap in Fig. [3] The dynamical instability at Q c i for low filling (n < 4/3) is therefore 
preempted by the closing of the single-particle excitation gap. 

We now turn to discuss the reason for the emergence of the flow-induced charge modulations in the kagomc lattice. 
To make the argument simpler, we restrict ourselves within the strong coupling regime (U 3> t) where fermion pairs 
become tightly bound molecular bosons. In this regime, the hopping term of the effective Hamiltonian for bosons 
in the presence of superflow Q is given by — JJ2(i j)( e ~ Z< ^ ' Tii b\t>j + H.c). Here, J = 2t 2 /U, b~i = e lCl ' ri c^cVf- is an 
annihilation operator for a boson, and = — rj is the bond vector. The Hamiltonian shows that the imposed 
superflow reduces the effective hopping amplitude of bosons by a factor of cos Q • ry , which plays a crucial role for 
the spatial modulation. The kagome lattice has three kinds of nearest-neighbor bonds connecting two sites: the A-B, 
B-C, and C-A bonds. The flow breaks the symmetry of the three bonds and introduces the anisotropy in the hopping 
amplitude, which naturally leads to the spatial modulation in density and order parameter. For example, when the 
flow is in the r — > K direction (Q = (Q, 0)), the effective hopping amplitudes are given by Jab = JcosQ for the A-B 
bond and Jbc = Jca = Jcos^ for the B-C and C-A bonds as shown in Fig. SJa). Because of this anisotropy, the 
system prefers forming a stripe modulation with ua = ^= iiq and Aa = Ab 7^ Ac in order to maximize the energy 
gain. Figure Etc) shows A y and n v as functions of the angle 6 of the flow momentum Q relative to the x axis. In 
general, A u and n v take different values depending on sublattices A, B, and C due to the anisotropy of the effective 
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FIG. 4: (color online) Effective hopping amplitudes for fermion pairs in the strong coupling regime when (a) Q = (Q,0) and 
(b) Q — (Q cos 9, Q sin 6). (c) Order parameter A^ and fermion filling n„ as functions of the angle 6 defined in (b). We set 
Q = 0.8tt, n = 1.6, and U/t = 3. 

0.5 
0.4 

o 

0.1 
0.0 




0.0 



0.5 1.0 1.5 

filling n 



FIG. 5: (color online) (a) Mean-field transition temperature T® as a function of fermion filling n. (b) Temperature dependence 
of order parameter modulation A, density modulation n, and averaged order parameter A. We set n = 1.6, Q = (0.6-7T, 0) and 
U/t — 3. The vertical dotted line marks T c ° in the absence of flow. 



hopping amplitudes of pairs, except for some special symmetric points, 9 = 0,7r/6,7r/3, • • • , where two of them are 
equivalent and the system forms a stripe pattern. This spatial modulation obviously emerges from the characteristic 
geometry of the kagome lattice and therefore it is absent in other typ ical lattice geometries such as square 14, lfj| . 
cubic jl4J,|l5|, honeycomb 17|,[27|, and even triangular lattice 2|f|29( that has geometric frustration. 

For realizing supcrfluid state in atomic fermions in optical lattices, experimental difficulties arise in cooling the 
system down to the supcrfluid transition temperature. The kagome lattice has a great advantage in this respect. 
Figure [5] (a) shows the mean-field transition temperature T c ° without flow as a function of the femion filling n. T c ° 
is significantly enhanced for high filling n > 4/3 when the Fermi level reaches the fiat band [3fjj . The experiments 
in Ref. [3 were performed in two dimension (2D). It may be necessary to cool the system below T c ° down to the 
Kosterlitz-Thouless transition temperature Ikt to achieve the superfluid quasi long range order in 2D. However, 
it is possible to suppress the phase fluctuations using a 3D optical lattice of weakly connected layered 2D kagome 
lattices. In this setup, the transition temperature is given by the mean-field value in Fig. [5] (a). Figure El[b) plots 
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the temperature dependence of h, A, and A for a superfluid state with stripe charge order. It shows that the stripe 
charge order is observable up to high temperatures slightly below T c °. We note that the curve for T c ° in Fig. [5] shows a 
dip in the vicinity of n = 2/3. This is due to the Dirac points in the free fermion band structure. The system is in the 
normal semimetallic phase for small U less than the critical value U c 2.8t at which the quantum phase transition 
to the superfluid phase takes place. The same phenomenon occurs in a honeycomb lattice [171 l3ll. |32|. 

The flow-induced charge modulations which we uncovered can be realized in cold atom experiments if superfiow is 
imposed by a moving optical lattice 12, 19|, 2(J. The experiment in Ref. [3 overlays two triangular lattices to form 



the kagome lattice. Each triangular lattice is formed by three lasers at angles of 120 degrees with respect to each 
other. Superfiow can be induced in the V — K direction by detuning one of the lasers for each triangular lattice and 
moving both the triangular lattices with the same velocity [lH ■ Superfiow can be also imposed by dipole oscillations 
which can be induced by a sudden displacement of a confining harmonic potential (2l| . Superfiow around a vortex is 
also expected to induce charge modulations. Since infinitesimally small flow can induce charge modulations, superfiow 
around a single vortex induces charge modulations that extend over the whole system. The charge modulations around 
a vortex may be therefore easily observed. 

In summary, we have studied the s-wave superfluid state of atomic fermions on the kagome lattice inspired by the 
recent experimental realization of tunable optical kagome lattices in Ref. Il8l We first perform a mean field analysis of 
the superfluid state while allowing for imposed superfiow. We find that superfiow induces a novel charge modulated 
state due to the characteristic geometry of the kagome lattice. The superfluid and CDW orders coexist in this state 
analogous to "supersolid" state. We examine the superfiow instabilities and critical current by evaluating the free 
energy as a function of superfiow. The critical current for high filling is found to be quite enhanced due to the flat 
band in the free fermion band structure. The superfluid state with charge modulations sustains high temperatures 
close to the transition temperature T c ° which is also enhanced by the flat band for high filling and therefore accessible 
using the setup of recent cold-atom experiments. 

We acknowledge T. Ohkane for performing preliminary calculations in the early stage of this work. S. T. and D. 
Y. were supported by Grant-in-Aid for Scientific Research, No. 24740276 (S. T.) and No. 23840054 (D. Y.). 
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Supplementary Material for "Flow-Induced Spatial Modulation of Superfluid Atomic 

Fermions in an Optical Kagome Lattice" 

A. The explicit form of ftq(k) 
The explicit form of the 6x6 matrix /icj(k) in Eq. (1) of the main text is given by 



'Q 



The matrix Kq can be diagonalized by the Bogoliubov transformation in a standard manner [l|. When the system 
does not have a superflow, we can take Aa = Ab = Ac = A and n\ = tib = nc = n due to the symmetry of the 
lattice. In this case, we can diagonalize /iq analytically, and the gap and number equations are easily derived as 

- vE E oeT— 7kT tanh o W 



( 


i 2 


"A 


-A A 


-2tcos(k+ • ^) 







-2tcos(k+ • 3f) 





\ 




-A A 




/i + fn A 





2icos(k_ • 


Si) 

2 > 





2£cos(k_ • !f) 






—2t cos(k + 


•¥) 





-fx - ^n B 


-A B 




-2tcos(k+ • ^) 














2icos(k_ ■ ^l) 


-A* 




3 





2tcos(k_ ■ ^ 






— 2t cos(k + 


•f) 





-2icos(k + • Sf) 







-fi - 


-A c 




\ 







2icos(k_ ■ 





2icos(k_ • 


Bl2 \ 

2 I 


— A* 


/i + ^n c 


/ 



U N ^ ^ 2E Q - 0T (k) 



and 



respectively. Here, N — 3M is the number of total lattice sites, f3 = l/T is the inverse temperature, and £ T (k) = 



e T (k) — /i — Un/2. The Bogoliubov quasiparticlc bands for Q = arc simply given by Eq =0 T (k) = y/£ T (k) 2 + |A| 2 , 
which means that there remains one fiat band (t = 0) even in the superfluid state (see Fig. 3 of the main text). 

B. Moving kagome optical lattice 

We present here how to prepare a moving optical lattice with kagome geometry. In the recent experiment by Jo 
et at, the kagome lattice was formed by overlaying two triangular optical lattices with different lattice constants 0]. 
Therefore we only have to consider the setup for moving a triangular-lattice potential with a constant velocity. The 
triangular lattice is generated by superposing three laser beams that intersect in the x-y plane with wave vectors kx 
= fc(l,0), k 2 = fc(— 1/2, — v3/2), and k3 = k(— 1/2, a/3/2). All beams are linearly polarized orthogonal to the plane 
and have the same field strength Eq. The total electric field is given by 

3 

E(r, t) = Y^Eo cos(k 4 • r - ut + <Pi)e z . (4) 

i=l 

The relative phases (f>ij = 4>i — <pj are fixed in Ref. to obtain a stable optical lattice. The generated dipole potential 
is proportional to the squared amplitude of the electric field 



|E tot | 2 - El 



3 + 2cos(bi • r + 23 ) + 2cos(b 2 • r + 3 i) + 2cos(b 3 • r + 4>i 2 ) 



2 

+ cos(2k! ■ r - 2ust + 2<t>i) + cos(2k 2 • r - 2ut + 20 2 ) + cos(2k 3 • r - 2ujt + 2cf> 3 ) 
+2cos((k 2 + k 3 ) • r - 2ut + (j> 2 + <f> 3 ) + 2 cos((k 3 + ki) • r - 2wt + </> 3 + </>i ) 
+2cos((ki+k 2 )-r-2o;i + 0i+0 2 )j, (5) 
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where = ei J 7 C (k J - — k/j). Since the frequency of light is quite large, only the time-averaged value of |E to t| 2 can affect 
atoms. Therefore, by dropping the terms containing 2ujt, we obtain a periodic dipolc potential 

V(r) = Vo Q + cos(bi • r + 2 a) + cos(b 2 ■ r + <j) 31 ) + cos(b 3 ■ r + 0i 2 )^ . (6) 

Red-detuned lasers give Vo < and we obtain a regular triangular-lattice potential, while the maxima of V(r) form 
a triangular lattice in the case of blue-detuned lasers (Vb > 0) 0]. Therefore, we can cancel out unwanted sites of a 
triangular lattice with Vo < by overlaying another potential with Vb > so that the total potential minima form a 
kagome lattice 0. 

One can move the lattice potential by introducing time-dependent phase differences </>y (t) through a small frequency 
detuning 5oj. Let us say that we detune one of the three lasers making a triangular lattice as (f>\ = Scot and <p2 = ^3 = 0. 
In this case, we can rewrite Eq. (|6]) as 

V(r) = Vb f ^ + cos(bi • r) + cos(b2 • r — Sut) + cos(b3 • r + Scot) 

= V ^ + cos (vW) + 2 cos Qfc (x + ^-tj ) cos ^j-ky] ^ , (7) 

which means that the potential moves in the T-K direction with a constant velocity 25uj/3k. A moving kagome optical 
lattice can be obtained by moving both the overlaid triangular lattices with a same velocity. Note that one uses two 
sets of lasers with different values of fc to create a kagome lattice Q . Therefore, the frequency detuning 6u> of each 
triangular-lattice potential has to be tuned so that the velocities 25u>/3k take a same value. 



[1] S. Tsuchiya, R. Ganesh, and A. Paramekanti, e-print [arXiv: 1206.4064vl. 

[2] G.-B. Jo, J. Guzman, C. K. Thomas, P. Hosur, A. Vishwanath, and D. M. Stamper-Kurn, Phys. Rev. Lett. 108, 045305 
(2012). 

[3] K. L. Lee, B. Gremaud, R. Han, B.-G. Englert, and C. Miniatura, Phys. Rev. A 80, 043411 (2009). 



